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Abstract 



The effective potential in the Grumiller's modified gravity [D. Grumiller, 
Phys. Rev. Lett 105, 211303 (2010)] includes the Newtonian potential and 
a Rindler term. The fitting to the rotation curve data of eight galaxies suggests 
a universal Rindler acceleration a « 0.30 x 10~ 10 m s -2 . We do a two-parameter 
fit first, with the mass-to-light ratio (T*) and the Rindler acceleration (a) as free 
parameters. It is found that the data of six out of the eight galaxies fit well with 
the prediction of theory in the range < r < 40 kpc, although the theoretical 
curves show a tendency of arising beyond this range. The Rindler accelerations 
of the six well-fitted galaxies have the same magnitude, with an average value 
a ~ 0.30 x 10~ 10 m s -2 . Inspired by this fact, we then carry out a one-parameter 
(T*) fit to the six galaxies, with a fixed at a, and find that the theory can still 
reproduce the observation. The value of Rindler acceleration we get here is a 
quarter of that of Milgrom's MOND. For the rest two galaxies, NGC5055 and 
DD0154, the fitting results are significantly improved if the photometric scale 
length (h) is included as another free parameter. 
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Introduction 



The observations on rotation velocities of stars around a galaxy center show significant 
discrepancies from Newtonian theory. According to Newtonian gravity, the rotation velocity 
is inversely proportional to the square root of distance from the center of a galaxy. However, 
the observed data often show an asymptotically flat rotation curve out to the furthest data 
points. Generally, there are three different ways to solve this problem. The most direct 
assumption is that there are a large number of nonluminous matters tha t have not been 
measured yet ( iBegeman et al.lll99ll ; iPersic et al.lll996t IChemin et l201lh . However, after 
decades of years heavy research, no direct evidences of the existence of dark matter have 
been found. This inspires the astronomers and physicists to search for other explanations of 
the discrepancy between the Newtonian dynamical mass and the luminous m ass. The most 
popular way is to modified the Newton dynamics (MOND) flMilgroml ll983al Jbl) . With only 
one universal parameter, ie., the critical acce leration q n , MOND has achieved great su ccess 
in explaining the mass discrepancy problem (jSandersI Il996l ; ISanders &: Verheijenl Il998l ) . 



Besides MOND, it is also possible to modify the Newtonian gravity (MOG). According 
to MOG, Newtonian gravity is inval i d at g alaxy scales. Such theories include n onsym 
metric gravity theory (NGT , Moffatl ( 119951 )). metric- skew-tensor gravity (MSTG , iMoffat 



(120051 ): iBrownstein fc Moffatl (J2006|)) 



Carmel ian general rela tivit y (CGR, 



, scalar- t ensor- vector gravity (STVG 



Moffat 



hood )) 



Carmelil ffioO^)), Horava-Lifshitz gravity theory (H-L 



theory, iHoraval J2009al lbllcl): ICardone etldl (bold ElJ)), etc.. All of them can in a large 



degree explain the observed rotation curves. 



Two years ago, iGrumillerl (120 lOf ) proposed an effective model for gravity at large dis- 
tance, and obtained an effective p otential in leading order c onsists of a Newtonian potential 
(oc 1/r) and a Rindler term (oc r). iGrumiller fc Preid (120111) showed that the Ri ndler force is 
capable of explaining about 10% of the Pioneer anomaly (lAnderson et al.lll998l ). and simul- 
taneously ameliorates the shape of galactic rotation curves. Most importantly, the Rindler 
term doesn't spoil the solar-system precision tests. 

In t his paper, we i nvestigate the galaxy rotation curves in the Grumiller's modified 
gravity (jGrumiller 2010 ). The rest of the paper is organized as follows: In section [21 we 
review briefly the Grumiller's modified gravity theory. In section [3J we do a best-fit procedure 
to constrain the parameters. The two-parameter fit shows that six out of eight galaxies have 
approximately the same Rindler acceleration a « 0.30 x 10~ 10 m s -2 . Then we carry out 
a one-parameter fit, with a fixed at a. It is found that the theory can still reproduce the 
observation. For the rest two galaxies, the fitting results are significantly improved if the 
photometric scale length is included as a free parameter. Finally, discussions and conclusions 
are given in section HI 
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2. Grumiller's modified gravity 



In Grumiller's modified gravity theory (iGrumillerl I201CH ) . the spacetime is described by 
a spherical symmetric metric 

ds 2 = g af3 dx a dx p + § 2 (d6 2 + sin 2 fl# 2 ), (1) 

where the two dimensional metric g a p (x 7 ) and the surface radi us <fe(x"<) depends o nly on the 
coordinates x 1 = {t,r}. The process of "spherical reduction" (IBerger et al.lll972l ) simplifies 
the 4-dimensional Eins tein-Hilbert action to a specific 2-dimensional dilation gravity model 
(IGrumiller et al.l 120021 ). The requirement of being power-counting renormalizable gives the 
action as 



K 



S = - d 2 x^g-[f($)R + 2(<9$) 2 - 2V($)\, 



(2) 



where k is a coupling constant, / and V are analytic functions of $. In order to go back to 
the Newtonian potential under some special conditions, we get / = $ 2 . 

To avoid curvature singularities for large $, V($) can't contain powers higher than $ 2 . 
The generic form of V($) reads as 



V{$) = A$ 2 + 5$ + b + £>(!/$). 



(3) 



By rescaling simultaneously $ and k, we can fix one of the coefficients. Without losing of 
generality, we choose 6 = 2. Substituting Eq.Q into Eq.(j2]) and dropping all subleading 
terms, and choosing convenient normalization of the coupling constants as well as k = 1, we 

get 

S = - 



d 2 x^f=g[<& 2 R + 2(<9$) 2 - 6A$ 2 + 8a$ + 2]. 
Differentiating the action in Eq.(jH) gives the equations of motion (EOM) 

R 



(4) 



4a 



ig af) V a d^ + 6A - 
2$(V„$, - g^V a d a )<S> - g^(d<S>) 2 = g^V^). 



(5) 



The general solution of the above EOM is 



' g afi dx a dx fi = -K 2 dt 2 + 
$ = r 



(6) 



where 



K 1 



2M , „ 
1 Ar 2 



2ar. 



(7) 
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If a = A = 0, Eq.(jSl) reduces to the Schwarzschild solution. While if M = A = 0, Eq.(j6]) 
reduces to the 2-dimensional Rindler metric. 

Consider a point particle with energy E and angular momentum I moving along a 
geodesic in the plane 9 = it/2. The conservation of energy and angular momentum gives 



r 1 + V cS = E, 
where 

T.eff M 12 Ml3 A ^ M l \ 

V = + W - — - T + Gr(1 + ^ (9) 
is the effective potential. The first to fourth terms on the right-hand- side of Eq.((9]) are the 
Newtonian potential, centrifugal barrier, general relativity correction and cosmic acceleration 
term, respectively. The last term, which is proportional to the Rindler acceleration, is 
novel. In the case of vanishing cosmological constant and vanishing angular momentum (ie., 
A = / = 0), the effective force corresponding to Eq.([9]) is given as 



rcff 



dV 



eff 



dr 



M 



a. (10) 

Z=A=0 r 



Thus, the rotation velocity of a testing particle moving in the potential of a point particle 
of mass M reads 



where v jv = ^M/r is the rotation velocity derived from Newtonian dynamics. This velocity 
has the asymptotic behavior 

v(r) = { (12) 
'or r — y oo. 



3. Best-fit procedure 

To calculate the theoretical rotation velocity, we should know the mass distribution 
of a galaxy. Generally, the mass in a galaxy contains two components: gas (mainly HI 
and He) and stellar disk. Some early type galaxies may also contain a bulge in the center. 
However, the mass of the bulge is often much smaller than that of the stellar disk, and we 
will neglect it for simplicity. Since the THINGS (The HI Nearby Galaxy Survey )0 brought 



1 http://ww.mpiarhd.mpg.de/THINGS/Data.html 
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an unprecedented level of precision to the measurement of the rotation curves of galaxies, 
we choose the THINGS galaxies as the samples. The profile of neutral hydrogen (HI) is read 
directly from the THINGS data cube (robust weighted moment-0) using the task ellint of 
Groningen Image Processing System (GIPSY]§ Assuming that HI locates in an infinitely 
thin disk, we calculate the rotation velocity contributed by HI (vm) directly using the task 

ROTMOD of GIPSY. 

The profile of stellar disk is derived from the photomet ric data. The surface brightness 
of the stellar disk can be fitted well by an exponential law (jde Vaucouleurslll959l ) 

/(r) = I(0)e-*, (13) 

where h is the scale length and 1(0) is the surface brightness at the center of galaxy. The 
integration of Eq. (fl3|) gives the total luminosity 

L = 2nl(0)h 2 . (14) 

Assuming that the mass-to-light ratio (Y*) is a constant and the disk is infinitely thin, the 
mass profile of the stellar disk can be written as 

S(r) = £(0)e - *, (15) 

where S(0) = Y*/(0) is the central surface mass density. The total mass of the stellar disk 
is given by 



M = 

The mass-to-light ratio is given as 



poo 

/ 2vrrS(r)rfr = 2ttS(0)/i 2 . (16) 
Jo 



t =M = m (17) 

* L 1(0)- [U) 
The total luminosity L of a galaxy is related to its absolute magnitude M. as 

M-A4 = -2.51og lo -^, (18) 

where A4 & and L Q are the absolute magnitude and total luminosity of the sum, respectively. 



The rotation velocity contributed by the exponential disk can be written as ( jFreeman 



197(1 ) 



M 

v,{r) = \l —l{r), (19) 



2 http://www. astro. rug.nl/ gipsy/ 
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where 



r 



3 



Here, 7 n and -fT n (n = 1, 2) are the n-th order modified Bessel functions of the first and 
second kind, respectively. 

The Newtonian velocity due to the combined contributions of gas and stellar disk is 
given by 

vn = \J\v 2 m + vl (21) 

where the factor 4/3 comes from the contribution of both helium (He) and neutral hydrogen 
(HI). Here we assume that the mass ratio of He and HI is Mne/Mni = 1/3. Any other 
gases are negligible compared to HI and He. Combining Eq.f llip and Eq. fl21l) . we get the 
theoretical rotation velocity in the Grumiller's modified gravity 



v(r) = y -v^(r) +vl{r) + ar. (22) 

In principle, Eq. (l22|) is only valid to point particle systems. Since the potential in Eq.([9]) 
is nonlinear to mass, superposition principle as well as Gauss theorem no longer holds, and 
we can't do an integration over the volume of mass as is done in the case of Newtonian theory 
to obtain the potential of any mass distribution. Any modified gravity are confronted with 
this problem. We adopt an approximation as is done in MOND and other modified gravities, 
by simply extrapolating Eq. (T22]) to any mass distribution, at least to the axisymmetric cases. 



Now we can use Eq. ( 122]) to fit the observed rotation curves. The free parameters are So 
(or Y*) and a. The sample galaxies and their properties are listed in Tabled] The rotation 
velocity data are read directly from the THINGS data cube (robust weighted moment-1) 
using the task rotcur of GIPSY, with PA, INCL, V sys and galaxy center (RA and DEC) 
fixed at the values given in Table [T] In order to constrain the parameters, we define the 
Chi-square as 

^W-fMF (23) 

i=l ai 

where w° bs is the observed rotation velocity, t> th (rj) is the theoretical velocity at radius 
calculated from Eq.( l22]) . and <jj is the uncertainty of v° hs . Then we employ the least-square 
method to minimize Eq. (1231) . The best-fit results are presented in Figure [1] and the values 
of parameters are listed in Table [2] 

From Table [2] we can see that the Rindler acceleration a has the same magnitude 
for all the eight galaxies except for NGC5055 and DD0154, with an average value a w 
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0.30 x 10~ 10 m s -2 . Inspired by this fact, we carry out a one-parameter fit procedure, with a 
fixed at a and S(0) as the only free parameter. The best-fit results are shown in Figure [21 
and the values of parameters are listed in Table From Figure [2J we find that the theory 
can still fits the observations (except for NGC3198), although the theoretical curves show a 
tendency of arising beyond the data range. For the two poorly-fitted galaxies, NGC5055 and 
DD0154, we include the disk scale length h as a free parameter and do a three-parameter 
fit. The results are plotted in Figure [31 and the best-fit parameters are listed in Table HJ 
It turns out that the fitting results are significantly improved. However, the best-fit scale 
lengths far deviate from the photometric scale lengths. We don't know if it mean that the 
mass scale length differs from that of photometry. 



4. Discussion and conclusion 

In this paper, we have done a best fit procedure to the rotation curves of eight THINGS 
galaxies in the framework of Grumiller's modified gravity. We did a two-parameter fit first, 
with So and a as free parameters. For the eight sample galaxies, six of them fit well in the 
range < r < 40 kpc. We found that the Rindler accelerations of each galaxies share the 
same magnitude, with an average value a w 0.30 x 10~ 10 m s~ 2 . The Rindler acceleration 
we got here is a quarter of the critical acceleration of Milgrom's MOND (do ~ 1.2 x 10~ 10 
m s -2 ). Inspired by this fact, we assumed that a may be a universal constant, and did 
another one-parameter fit to the six well-fitted galaxies, with a fixed at a and Y* as the 
only free parameter. It turned out that the theoretical rotation curves can still reproduce 
the observed data, with only one exception. The theoretical velocity of NGC3198 shows a 
tendency of sharply arising after ~ 20 kpc. For the rest two galaxies, NGC5055 and DD0154, 
if we include the disk scale length as a free parameter, the fitting results are significantly 
improved. Thus, the Rindler acceleration of all the eight galaxies, except for DD0154, can 
be unified to a. However, the Rindler acceleration of DD0154, which is a gas rich dwarf 
galaxy, is one order of magnitude smaller. 

In t he fitting proced ure, we regarded the Rindler acceleration as a constant for a certain 



galaxy. iGrumillerl (120101 ) argued that a may varies from system to system. Further more, 
a is not necessary to be a constant even in the same system: it may be a function of r. 
For instance, a is one order of magnitude larger in the Sum-Pioneer system than that of in 
the Galaxy-Sum system. In the Earth-Satellite system, it may be much larger. Otherwise, 
the theory can't reconcile with all the experiments. Unfortunately, it is still unknown what 
determines the scale of a. In this paper, however, we found that a may be a universal 
constant. It is surprising that the value of a we obtained here is same to the upper bound 
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value constrained from the Cassini spacecraft data (IGrumiller fc Preisl l201ll ) 



The most interesting feature of Grumiller's gravity is that it predicts a rotation curve 
proportional to the square root of distance in the large distance limits. However, the observed 
data just extend to a few tens of kpc. Although the velocity law Eq. (l22[) follows well to the 
experiment data in this range, it's too early to say that this is a valid model at large distances. 
Future observations of rotation velocity at large distances may provide further tests to the 
theory. It is very likely that Gru miller's gravity w ill fail in the large distance limits if a is 
a constant. As was mentioned by iGrumillerl (120101 ). a may be a function of r. If we choose 
a = y/GMao/r, Eq. ffTTj) becomes 



+ y/GM{ 



r)a , 



(24) 



where ao ~ 1.2 x 10 10 m s 2 is the critical acceleration in Milgrom's MOND and M(r) is 



the mass surrounded by a sphere of radius r. Eq.(l24l) has the sa me asymptotic behavior as 
MOND, with the Tully-Fisher relation holds. iLi fc Chang obtained a similar formula 

with the same asymptotic behavior as Eq. (1241) based on the Finsler geometry. ILi et al.l (120121 ) 
showed that Eq.(!24l) plus the dipole and quadrupole contributions can well explain the mass 
discrepancy of Bullet Cluster 1E0657-558. Since the observed data show that most galaxies 
have asymptotically flat rotation curves, it may better describe the galaxy rotation curves if 
a is inversely proportional to r than is a constant. 



We are grateful to Y. G. Jiang and S. Wang for useful discussion. This work has been 
funded in part by the National Natural Science Fund of China under Grant No. 10875129 
and No. 11075166. 
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Table 1. Properties of the sample galaxies. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Names 


RA 


DEC 


INCL 


PA 


^sys 


D 


h 


M m 


m B 


M B 


L 




[h m s] 


[° ' 1 


[°] 


[°] 


[km/s] 


[Mpc] 


[kpc] 


[1O 8 M ] 


[mag] 


[mag] 


[1O 1O L ] 


NGC2403 


07 36 51 1 


+65 36 03 


63 


124 


132.8 


3.2 


2.05 


25.8 


8.11 


-19.43 


0.920 


NGC2841 


09 22 02.6 


+50 58 35 


74 


153 


633.7 


14.1 


3.55 


85.8 


9.54 


-21.21 


4.742 


NGC2903 


09 32 10.1 


+21 30 04 


65 


204 


555.6 


8.9 


2.81 


43.5 


8.82 


-20.93 


3.664 


NGC3198 


10 19 55.0 


+45 32 59 


72 


215 


660.7 


13.8 


3.88 


101.7 


9.95 


-20.75 


3.105 


NGC3521 


10 05 48.6 


-00 02 09 


73 


340 


803.5 


10.7 


2.86 


80.2 


9.21 


-20.94 


3.698 


NGC5055 


13 15 49.2 


+42 01 45 


59 


102 


496.8 


10.1 


5.00 


91.0 


8.90 


-21.12 


4.365 


NGC7331 


22 27 04.1 


+34 24 57 


76 


168 


818.3 


14.7 


4.48 


91.3 


9.17 


-21.67 


7.244 


DD0154 


12 54 05.9 


+27 09 10 


66 


230 


375.9 


4.3 


0.54 


3.58 


13.94 


-14.23 


0.008 



Note. — Column (1): galaxy names. Columns (2) and (3): galaxy centers in J2000.0 from I Walter et al.l (120081). Columns 
(4) an d (5): Inclinations and positio n angles from IWalter et al.1 (120081 ). Column (6): systematic velocities from de Blo k et al 



2008 ) . Co lumn (7): dis tances from I Walter et al.l ( 20081 ) . Column (8) : the scale lengths of optical disk from Begeman et al 



1991 ) and iFlored (|1993l ) , but corrected with the new dis tances in column (7 ) . Column (9) : mass of HI from 



Walter et al 



2008) . Column ( 1 0): a pparent B-band magnitudes from IWalter et al.l ((2008). Column (11): absolute B-band magnitudes 



from lWalter et al.l (|2008l ). Column (12): luminosity calculated from column (11) using Eq. (|18|) . 
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Table 2: The two-parameter best fit. 





So 


a 




M/L 


X 2 /n 




[M & pc- 2 ] 


[10- 10 m s- 2 ] 


[1O 1O M ] 


[M e /L Q ] 




NGC2403 


450.75 ± 9.73 


0.301 ±0.006 


1.19 ±0.03 


1.29 ±0.03 


1.21 


NGC2841 


2626.12 ±269.17 


0.396 ±0.056 


20.79 ±2.13 


4.39 ±0.50 


2.47 


NGC2903 


1190.88 ±24.66 


0.316 ±0.010 


5.91 ±0.12 


1.61 ±0.03 


3.87 


NGC3198 


421.25 ±9.75 


0.160 ±0.005 


3.98 ±0.09 


1.28 ±0.03 


2.49 


NGC3521 


1239.17 ±40.85 


0.296 ±0.014 


6.37 ±0.21 


1.72 ±0.06 


0.54 


NGC5055 


900.27 ± 14.72 


0.053 ±0.013 


14.14 ±0.23 


3.24 ±0.05 


5.84 


NGC7331 


1049.41 ± 53.30 


0.299 ±0.021 


13.23 ±0.67 


1.83 ±0.09 


0.35 


DD0154 


92.76 ±46.23 


0.060 ±0.010 


0.02 ±0.01 


2.12 ± 1.06 


17.48 



Note. — The free parameters are So and a. M^^ is the mass of the disk calculated from Eq. (|16p . 
M/L is the mass-to-light ratio calculated from Eq. (fT7|) . x 2 / n is the reduced chi-square. The 
numbers after "±" are the la errors of the corresponding parameters. 
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Table 3: The one-parameter best fit. 





So 


M disk 


M/L 


X 2 /n 




[M pc- 2 ] 


[1O 1O M ] 


[M /L ] 




NGC2403 


444.96 ± 10.78 


1.17±0.03 


1.28 ±0.03 


1.23 


NGC2841 


2825.29 ±93.10 


22.37 ±0.74 


4.72 ±0.16 


4.15 


NGC2903 


1209.00 ±34.15 


6.00 ±0.17 


1.64 ±0.05 


3.89 


NGC3198 


248.29 ± 13.20 


2.35 ±0.12 


0.76 ±0.04 


18.79 


NGC3521 


1222.09 ± 64.39 


6.28 ±0.33 


1.70 ±0.09 


0.55 


NGC5055 


X 


X 


X 


X 


NGC7331 


1037.32 ± 72.80 


13.08 ±0.92 


1.81 ±0.13 


0.36 


DD0154 


X 


X 


X 


X 



Note. — The free parameter is T,q. The parameter a is fixed at the average value a ~ 0.30 x 
10~ 10 m s -2 . is the mass of the disk calculated from Eq. (|16p . M/L is the mass-to-light ratio 

calculated from Eq. (jl7j) . x 2 /n is the reduced chi-square. The numbers after "±" are the la errors 
of the corresponding parameters. We skip the one-parameter fit for NGC5055 and DD0154. 



Table 4: The three-parameter best fit. 





So 


a 


h 


M disk 


M/L 


X 2 /n 




[M e pc- 2 ] 


[10- 10 m s- 2 ] 


[kpc] 


[1O 1O M ] 


[M /L ] 




NGC5055 


1503.24 ±217.05 


0.245 ± 0.042 


2.81 ±0.32 


7.46 ± 2.78 


1.71 ±0.64 


2.50 


DD0154 


47.86 ±3.87 


0.027 ±0.005 


2.45 ± 0.28 


0.18 ±0.06 


22.56 ±6.98 


0.86 



Note. — The free parameters are So, a and h. M^sk is the mass of the disk calculated from 
Eq. (|16p . M/L is the mass-to- light ratio calculated from Eq. (|17p . x 2 / n is the reduced chi-square. 
The numbers after "±" are the la errors of the corresponding parameters. 
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NGC2403 NGC2841 




Fig. 1. — Two-parameter fit (black curves) to the rotation curves of sample galaxies. The x-axis 
is the distance in kpc, and the y-axis is the rotation velocity in km s _1 . The red curves are the 
contribution of gas (HI and He). The green curves are the contribution of gas and stellar disk in 
Newtonian dynamics. 
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Fig. 2. — One-parameter fit with a fixed at a ~ 0.30 x 10~ 10 m s~ 2 . The x-axis is the distance in 
kpc, and the y-axis is the rotation velocity in km s _1 . The meaning of each line is the same to that 
in Figure [TJ 

NGC5055 DD0154 




5 10 15 20 25 30 5 10 15 



Fig. 3. — Three-parameter fit with So, a and h free. The x-axis is the distance in kpc, and the 
y-axis is the rotation velocity in km s -1 . The meaning of each line is the same to that in Figure [TJ 



